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ABSTRACT 

We developed a new model for the X-ray spectral fitting xspec package which takes into account 
the effects of both thermal and dynamical (i.e. bulk) Comptonization. The model consists of two 
components: one is the direct blackbody-like emission due to seed photons which are not subjected to 
effective Compton scattering, while the other one is a convolution of the Green's function of the energy 
operator with a blackbody-like seed photon spectrum. When combined thermal and bulk effects are 
considered, the analytic form of the Green's function may be obtained as a solution of the diffusion 
Comptonization equation. Using data from the BeppoSAX, INTEGRAL and RXTE satellites, we 
test our model on the spectra of a sample of six persistently low magnetic field bright neutron star 
Low Mass X-ray Binaries, covering three different spectral states. Particular attention is given to the 
transient powerlaw-like hard X-ray 30keV) tails that we interpret in the framework of the bulk 
motion Comptonization process. We show that the values of the best-fit (5-parameter, which represents 
the importance of bulk with respect to thermal Comptonization, can be physically meaningful and 
can at least qualitatively describe the physical conditions of the environment in the innermost part 
of the system. Moreover, we show that in fitting the thermal Comptonization spectra to the X-ray 
spectra of these systems, the best-fit parameters of our model are in excellent agreement with those 
of COMPTT, a broadly used and well estabhshed xspec model. 

Subject headings: stars: individual: Sco X-1, GX 17-1-2, Cyg X-2, GX 340-hO, GX 3+1, GX 354-0 — 
stars: neutron — X-rays: binaries — accretion, accretion disks 



1. INTRODUCTION 

The study of the transient hard X-ray tails in neu- 
tron star (NS) Low Mass X-ray Binaries (LMXBs) has 
recently received a strong theoretical impulse. The first 
approach has been mainly phenomenological, with the 
hard tails simply fitted by a powcrlaw (pl) model (see 
review by Di Salvo & Stella 2002). Subsequently, using 
0.4-120 keV BeppoSAX ddXa., FarineUi et al. (2005, here- 
after F05) tried a more physical approach, describing the 
hard tail of GX 17+2 by means of Comptonization of soft 
(~ 0.6 keV) seed photons off a hybrid (thermal plus non- 
thermal) electron population. More recently, Paizis et 
al. (2006, hereafter P06), using long-term time-averaged 
INTEGRAL/lBlS spectra (> 20keV), confirmed the 
presence of hard X-ray emission in the Z sources Sco X-1, 
GX 5-1, GX 17+2, GX 340+0, Cyg X-2, and they also 
discovered it in the atoll source GX 13+1, which, in fact, 
is characterized by a stable Z-like X-ray continuum (ther- 
mal component with low electron temperature kT^ and 
high optical depth t). 

Motivated by observations of the intermediate and 
soft states of black hole candidates (BHCs), P06 pro- 
posed that bulk motion Comptonization can be respon- 
sible for the formation of the well-known hard X-ray 
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tails in NSs observed in these sources as well. In their 
analysis, they used the Bulk Motion Comptonization 
(bmc) model in xspec, whose description is reported 
in Titarchuk et al. (1996, 1997, hereafter TMK96 and 
TMK97, respectively). Moreover, merging their results 
with those of Falanga et al. (2006, hereafter F06) on the 
atoU source 4U 1728-34 (GX 354-0), P06 identified four 
main spectral states for NS LMXBs (see Fig. 4 in P06): 
hard/PL state (H 1750-440 and H 1608-55), low/hard 
state (GX 354-0), intermediate state (e.g., GX 5-1) and 
very soft state (e.g., GX 3+1). In their picture, the rel- 
ative contribution of thermal and bulk Comptonization 
(TC and BC, respectively), in turn related to the accre- 
tion rate M and the radiation pressure from the NS, is 
the main parameter explaining both the different spec- 
tral states among several sources and the spectral evolu- 
tion of a single source. However, the lack of data below 
20keV prevented P06 from drawing more stringent con- 
clusions on the accretion geometry of the systems and on 
the physical parameters. 

Re-analyzing broad-band BeppoSAX data, Farinelli et 
al. (2007, hereafter F07) appfied the BMC model to 
GX 17+2. They found that the source spectrum can be 
described by a blackbody (BB)-like emission, plus an un- 
saturated TC component plus a PL. The TC component, 
fitted with the COMPTT model (Titarchuk 1994, hereafter 
T94) , is suggested to come from the region (called transi- 
tion layer, TL) between the Keplerian accretion disk and 
the NS surface, while the PL-like emission is attributed 
to BC of the NS seed photons on the in-falling mate- 
rial, which is in turn the innermost part of the TL itself. 
This accretion scenario proposed by F07 is supported 
by results obtained by numerical solutions of the equa- 
tion for the radial velocity profile in accretion disks 
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(Titarchuk & Farinelli 2008, hereafter TF08, in prepara- 
tion), where the quasi free-fall behaviour of wr inside the 
TL from some radius Rg to the NS surface, is unambigu- 
ously shown. 

It is thus evident that BC is gaining strong theoretical 
and observational support in the study of LMXB spec- 
tral evolution. This is true also for the case of accretion 
powered X-ray pulsars for which a new theoretical model 
based on TC and BC, occurring in the accreting shocked 
gas, was recently presented (Becker & Wolff, 2007, here- 
after BW07). In the BW07 Comptonization model the 
effects of the strong magnetic field {B ^ 10^^ G) are very 
important, as the accreting gas column is channelized at 
the NS magnetic poles and the presence of this field also 
has important consequences for the photons propagating 
through the plasma. 

From the observational point of view, the energy index 
a and the cutoff energy are only parameters which can 
be determined by fitting the observable spectral shapes 
to the Comptonization models; additional assumptions 
have to be adopted however in order to relate them to 
the physical parameters of the particular models. 

However in this Paper we concentrate in the study of 
sources where the magnetic field is low, so that it has 
second-order effects with respect to radiation pressure 
in driving the accretion process. For a such class of 
low-magnetized compact objects, in the current xspec 
package it is available the BMC model, which was used, 
as mentioned above, by P06. It is in fact a general 
Comptonization model, i.e. widely applicable, given that 
the output emerging Comptonized spectrum is given by 
the convolution of the analytical approximation of the 
Green's function (response of the system to injection of 
a monochromatic line at energy Eq) with a BB seed pho- 
ton spectrum, regardless of the Comptonization process, 
i.e. only due to thermal or combined thermal and bulk 
effect. 

The BMC model parameters are the bb seed photons 
color temperature (/cTbb) and normalization, the energy 
index a of the Comptonized spectrum and the illumina- 
tion factor log(yl.). The information on the efficiency 
of the Comptonization is carried-out by the index a: 
the lower a, the higher the energy exchange from hotter 
electrons to softer seed photons. As we mention above, 
a does not specify which kind of Comptonization pro- 
cess produces it, but it is simply related to an observ- 
able quantity in the photon spectrum of the data. The 
Green's function of the model is approximated by a bro- 
ken PL, a condition that holds when the average energy of 
the seed photons is much lower than the average energy 
Eav of the plasma (Sunyaev & Titarchuk 1980, hereafter 
ST80, and TMK97). The BMC model however lacks the 
electron recoil-effect term in the Green's function, which 
is responsible for the rollover in the spectrum (the energy 
of the rollover, depending on the plasma temperature 
and on the first and second order of the bulk velocity, at 

about i?av)- 

Thus the Comptonized part of the emerging spectrum 
is in the BMC model always a PL with no cutoff. How- 
ever the lack of the energy cutoff in the model becomes 
critical in fitting of the TC spectra of LMXBs in their 
high-luminosity state, when the rollover is well observed 
around lOkeV. 

Motivated by these limitations in the currently avail- 



able version of the BMC model for xspec and by the 
recent observational results on LMXBs (see above), we 
developed a new model for XSPEC which can be consid- 
ered as an extension and completion of BMC. In Sec- 
tion [2] we summarize the main results of the Bulk Mo- 
tion Comptonization theory and show the mathematical 
background of our model. In Section [3] we show the re- 
sults of the X-ray spectral fitting using our model on a 
sample of six persistently bright NS LMXBs. In Section[4] 
we discuss the main results and the limits of applicability 
of the model itself, with a breif theoretical comparison 
to the model of BW07, and in Section [5] we draw our 
conclusions. 

2. DESCRIPTION OF THE MODEL 

2.1. Solution of the Comptonization equation for bulk 
flow of non-relativistic electrons 

The Comptonization equation (a particular type of 
Fokker-Planck equation, see e.g. T94) describes an evo- 
lution of the photon field due to Compton scattering 
of photons off non-relativistic thermal electrons. It is 
obtained as the approximation of the full kinetic Boltz- 
mann equation when the photon field is almost isotropic, 
the average number of scatterings is high and the av- 
erage energy exchange per scattering is small, namely 
AE/E <C 1. When the divergence of the velocity field 
V • Wb of the electrons is zero (namely when electrons 
are subjected to pure Brownian motion) only TC con- 
tributes to the total emerging spectrum (ST80). On the 
other hand, if the plasma flow is subjected to inward 
bulk motion (V • ?;b < 0), the latter provides a further 
channel to Comptonization. The Comptonization equa- 
tion including the bulk motion term was first derived 
and introduced by Blandford & Payne (1981). A slightly 
modified version of the equation, which includes the sec- 
ond order bulk term Vh^ is presented in TMK97. The 
equation for the photon occupation number n(x, t) in 
the case of a free-fall velocity profile {vr oc r~^/^) reads 
as 
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(1) 



where x = E/kT^ is a dimensionless energy, O = 
kTc/mcC^ is a dimensionless electron temperature, to 
is the accretion rate in units of the Eddington accre- 
tion rate, k is the inverse of scattering free path of 
electron and r is the optical depth. The term j in 
equation ([T]) is a function of x and r and represents 
the seed photons source function and S^^ = S^^fbir), 
/,(r) = l+(yc)7(3e). 

The quantity of interest in the above equation is the 
bulk parameter S, which is derived in TMK97 as 



Or 



(2) 



where i = L/LEdd is the fractional Eddington luminosity 
impinging on the flow. 
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It is worth noting that in equation ([T]), the term 
n depends on r through the bulk term v^. Given 
that ff^^n depends on r it is possible to find an an- 
alytical solution of the equation (using the separation 
of variables method) only neglecting the v'^/c^ term, 
namely by setting fb—l- In this case the Green's func- 
tion G{x,xo)=Ng{xo)G{x,xo) (see also Eqs. [B4]- 
in TMK97) with G{x,xo) given by 

G{x,xo) 



xa r(2c 



S) 



, xQ+5+3 

iFi(a,4 + 2a + (5,x) J(a;o, a, w), x < Xq; 

iFi(a,4 + 2a + (5,Xo) J(x,a,t^), x>Xq. 

Here a is energy index, ll'=q;-|-(5-|-3 and x, xq represent 
the dimensionless scattered and injected photon energies, 
respectively, namely x = E/kTc, xq = Eo/kTc. This 
is in fact the complete form of the Green's function for 
the Comptonization equation, while in equation (22) of 
TMK97 it is reported the asymptotic case of low-energy 
(a;o ^ 1) monochromatic line injection. The dependence 
of spectral energy index a on the model physical param- 
eters O, T and 6 is presented in equations (7) and (23) 
of TMK96 and TMK97, respectively. We note that the 
dependence of a on O and t in equation ([3]) is implicit 
through the presence of the first eigenvalue of the space 
diffusion operator. 

In equation iFi is a confluent hyper-geometric 
function (Abramowitz & Stegun 1970), and J(a;,Q;,aj) 
is an integral function expressed through the formula 



dt. 



(4) 



J{x,a,uj) = / e^'{x + t)'^ t° 
Jo 

Given that Comptonization conserves the photon num- 
ber, this photon conservation law can be rewritten as the 
integral equation for the Green's function: 



G{x,xo)^^^^ 



(5) 



/o a; Xq 
Thus the normalization factor Nq{xq) in equation ([3]), 
which is dependent of the injected line energy xq reads 
as 



Ng{xo) - 



Xo 



G{x,xq) 



dx 



(6) 



When the energy of the injected line is much smaller than 
the plasma energy, i.e. when xq, x uj, the Green's 
function is a broken PL and iVG(a;o) ~ a{a -I- 5 -t- 3) (see 
TMK97). 

The emergent Comptonized spectrum is given by the 
convolution of the Green's function with S{xo) if the in- 
put spectrum, instead of being a simple monochromatic 
line, has a continuum shape S{xq). Let us consider seed 
photons distributed according to the law 



S{xq) 



(7) 



where and Tg are the electron and seed photon tem- 
peratures, respectively, and Cn is the normalization fac- 
tor. This is evidently a blackbody spectrum if 7=3. The 
Comptonization spectrum is thus given by 




Fig. 1. — The Green's functions G{x,xo) for monochromatic hne 
with input dimensionless energy a;o=0.3 and different values of a 
and S: a=l, 5=0 (black); a=l, <5=30 (red); a=2.5, S=0 (green); 
o=2.5, (5=30 (blue). The increase of the cutoff energy due to bulk 
effect is clearly visible. 
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dxo. (8) 



Given that the two integrals in equation ^ do not have 
analytic presentation, we have to calculate them numer- 
ically. This step is performed using the Gnu Scientific 
Library^ package. Note also that in the same equation 
the functions Ng{xo) and J(xo,a,uj) are themselves in- 
tegrals (see Eqs. [3] and [6j ) and thus should be prelim- 
inary calculated before putting them into the integral. 
This embedded integral operation is however very high 
CPU-time consuming, and it makes the whole convolu- 
tion process very slow, for xspec running. 

To avoid this problem, the integral function in equa- 
tion ([4]) is analytically solved at each integration step us- 
ing the steepest descend method^ described in Appendix 
C of TMK97. This method provides very high accu- 
racy (within few percent) for high (5- values (> 10), while 
when S is close to zero the discrepancy may reach 10% 
(depending on the corresponding a and a;- values). To 
better approximate the integral value when S is low, we 
produced a three-dimensional correction matrix of S, a 
and x-values (with S in the interval 0-10) which provides 
a correction factor between the numerical and analytical 
results of integral ([4]). The correction factor is computed 
using linear interpolation among the 6, a and x-bounds 
of the grid. The same approach is adopted to compute 
the normalization constant Ng{xo); for (5 > 10 we simply 
use the value a{a + S + 3), while for lower (5-values we 
use an interpolation grid. A subsequent numerical check 

® The Gnu Scientific Library is freely available at the web site 
www.gsl.org 
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reveals that the photon number is conserved with high 
accuracy (within few percent) over a very large spanning 
of the a and ^- values. 

2.2. The free parameters of COMPTB 

Following the same approach of the BMC model 
(TMK96, TMK97), we structure our model as being the 
sum of two components: one represents the emission of 
the soft seed photons which are not affected by a notice- 
able up-scattering in the plasma cloud, while the other 
one is a result of efhcient Comptonization. The emergent 
spectrum, similarly to BMC, can thus be described as 



F{E) 



A + l 



(BB + Ax BB*G), 



(9) 



where BB = S{xo) (see Eq. [7]) and the convolution 
BB * G is given by equation ([S]). The free parameters 
of our model are the seed photon temperature kTs, the 
index 7 (if 7=3 seed photons have pure BB spectrum), 
the Green's function energy index a (see Eq. [3]), the 
bulk parameter 6, the temperature of the electrons kTc, 
the illuminating factor log(A) and the normalization con- 
stant Cn. The latter is chosen is such a way that when 
A — > (and 7=3) the model simply reduces to the stan- 
dard BB model for xspec. We note also that Cn must be 
multiplied by the factor 0.5 x 10'^^(kTc/'keV)^ in order to 
express the flux in physical units (ergcm~^ s~^ keV^^). 
In equation ((9|) the factor 1/(1 + A) is the fraction of 
the seed photon radiation directly seen by the observer, 
whereas the factor A/{1 + A) is the fraction of the seed 
photon radiation up-scattered by the Compton cloud. 

As it can be clearly seen, we expand the parameter 
set with respect to BMC. In principle, this allows one to 
extract a more detailed physical information from the 
observed spectra. In Figure [T] we demonstrate the bulk 
motion effect presenting the Green's functions G{x,xo) 
for a couple of a-values (1 and 2.5) and two different 
(S-values (0 and 30) for each a. The cutoff energy is 
higher for greater (5- values (see red and blue with respect 
to black and green curves related to G(x,xo) for which 
6=0). Note that the photon number Nc{xo) is the same 
for all G{x,Xq). We want also to emphasize however, 
that our model is not specific to bulk motion: it is evident 
indeed from equations Q and ([9|) that setting ^ ^ 1 and 
(5=0, a pure TC spectrum is reproduced, like the mostly 
used models compst (ST80) and comptt (T94). 

It is also important to point out the difference between 
the standard ones and our Comptonization model. In 
our case the code provides, as the best-fit parameters, 
the electron temperature kTc and the energy spectral 
index a instead of kT^ and optical depth r as in COMPST 
and COMPTT. Thus, once kTg and a are provided, it 
is possible to infer the Thomson optical depth r using 
the relations reported in Titarchuk & Lyubasrkij (1995, 
hereafter TL95) for both spherical and slab geometry 
of the Comptonizing plasma. For example, in the non- 
relativistic limit the following equation holds: 




Cr kTc (t + 2/3)2 



(10) 



where r is the optical radius of the sphere and optical 
half-thickness of the slab for spherical and plane geome- 
tries, respectively, Ct=3 for a sphere and Cr=12 for a 



slab. When bulk is present (5 > 0), the relation between 
a, S and kTc is not straightforward and must be calcu- 
lated numerically (see TMK97). 

3. APPLICATION OF THE MODEL TO A SAMPLE OF 
LMXBS 

In order to test our model, we used a data set from a 
sample of six LMXBs belonging to the Z and atoll classes 
(Hasinger & van der Klis 1989): sources of the former 
class are Sco X-1, GX 17+2, Cyg X-2 and GX 340+0, 
while GX 354-0 (4U 1728-34) and GX 3+1 belong to the 
latter one. 

The choice of the first four sources was mainly moti- 
vated by the fact that they show transient PL-like hard 
(> 30keV) X-ray emission over their stable continuum 
(Sco X-1: D'Amico et al. 2001, P06; GX 17+2: Di Salvo 
et al. 2000b, F05, F07; Cyg X-2: Di Salvo et al. 2002, 
hereafter DS02; GX 340+0: Lavagetto et al. 2004, here- 
after L04). Other sources, such as GX 5-1 (Asai et al. 
1994, Paizis et al. 2005), GX 349+2 (Di Salvo et al. 
2001) and, most recently, GX 13+1 (P06) have shown 
the same high-energy transient behaviour. The stable 
continuum spectra of these source are usually fitted with 
a two-component model consisting of a photoelectrically- 
absorbed BB plus the TC model COMPTT by T94; when 
the hard X-ray tail appears, an additional PL is included 
in the model to fit to the data. The presence of a gaus- 
sian emission line around 6.7 keV is also observed in these 
sources. 

The persistent continuum spectrum of the bursting 
atoll GX 3+1 is very similar to that of classical Z sources, 
with a Comptonization spectrum characterized by kTc ~ 
3keV and r - 10 (Oosterbroek et al. 2001, P06) plus 
contribution from a ~ 1 keV BB-like component. Nei- 
ther spectral transitions nor hard X-ray tails have been 
observed up to now in this source. 

Finally, the bursting atoll source GX 354-0 is char- 
acterized by the presence of different spectral states 
that correlate with its position on the hardness-intensity 
diagram (HID). Long-term monitoring performed with 
INTEGRAL by F06 has shown that in the lower island 
state of the HID, the source spectrum is quite hard and 
can be described by a TC spectrum with kTc ~ 30keV 
and optical depth r ~ 1.5 (assuming slab geometry) with 
no evidence of soft BB-like emission. As the source moves 
across the HID, reaching the banana state, the TC spec- 
trum gets softer, with electron temperature dropping to 
~ 3 keV (with r correspondingly increasing to ~ 5) and 
simultaneous appearance of a soft (~ 0.6 keV) BB-like 
component which was attributed by F06 to the accretion 
disk. 

Selection of a sample of sources like these thus al- 
lows us to test our model over a wide range of spec- 
tral states which are experienced by NS LMXBs, as it 
is also noticed by P06. Following the state terminol- 
ogy of P06, we can apply our model to the interme- 
diate state spectra (Sco X-1, GX 17+2, Cyg X~2 and 
GX 340+0) where the combined effect of the thermal 
and bulk Comptonization is suggested to be the origin 
of the hard X-ray emission, to very soft state spectra 
(GX 3+1), where thermal Comptonization by a rela- 
tively cold plasma {kTc ~ 3 keV) dominates, and to 
low/hard state spectra (GX 354—0), in which the thermal 
Comptonization component characterized by relatively 
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Fig. 2. — Schematic view of the proposed geometry for thermal and bulk Comptonization regions in LMXBs hosting a NS with PL-like 
emission at high energies. The thermal Comptonization spectrum (thermal comptb) originates in the outer part of the TL region, whereas 
the thermal plus bulk Comptonization spectrum (thermal plus bulk COMPTb) arises in the innermost part of the TL where the BB-like (NS) 
seed photons are (thermally and dynamically) Comptonized by the in-falling material. The radial extension of the two Comptonization 
regions is, in fact, rescaled for clarity, whereas numerical calculations show that they are of order -Rns- 



hot electron temperature {kTc ~ 30keV) dominates. 

We refer the reader to P06 for a more complete discus- 
sion on the interpretation of the evolution among these 
spectral states. It is worth noting that the fit of the 
fourth spectral state identified by P06 {hard/PL state) 
with COMPTB (similarly to COMPTt) will not provide 
physically meaningful parameters, due to the presence of 
the non-attenuated PL that results in very high plasma 
temperatures (hundreds of keV), i.e. out of the regime of 
validity of the model. The spectral shape can be, in fact, 
fitted, but the resulting parameters are not physically 
meaningful (see discussion in Section 4.3). 

As a complement to the interpretation by P06, in their 
analysis of GX 17+2, F07 proposed an accretion scenario 
in which the BB-like component in the X-ray spectrum 
arises from the NS surface, while the TC spectrum orig- 
inates in the outer part of the TL region. The PL-like 
hard X-ray emission is instead the result of bulk Comp- 
tonization of the BB-like (NS) seed photons by the inward 
in-falling material in the innermost part of the TL (see 
Fig.©. 



In the analysis of the intermediate state sources, we 
proceed in the following way: first we fit their spectra 
with a mostly used multi-component model consisting 
of a BB plus COMPTT plus a PL plus, when required, 
GAUSSIAN emission lines (see above). Given that the ra- 
dial extension of the TL is small compared to its ver- 
tical height-scale (e.g., Titarchuk & Osherovich 1999, 
TF08), so that it can be somewhat approximated to a 
geometrically thick equatorial belt around the NS, we 
assumed a slab geometry for COMPTT for all the sources 
but GX 354-0, whose spectral state seems to be more 
related to a quasi-spherical accretion. Using this multi- 
component model we do not report the unabsorbed bolo- 
metric (0.1-200 keV) source fluxes because the presence 
of the PL component gives origin to a strong bias at low 
energies, where in fact the flux diverges, and this effect 
gets worse as the PL steepness increases. It is worth 
pointing out this problem when one fits spectra with PL 
components. 

Subsequently, we replace COMPTT with COMPTB, fix- 
ing the illuminating factor to 1 (i.e. log(74)=8, cor- 
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SCO X-1 GX 1 7+2 




Energy (keV) Energy (keV) 

Fig. 3. — Unabsorbed EF(E) spectra of the six studied NS LMXBs, superposed with the best-fit models (solid line, see Table [2]l and 
residuals between the data and the models in units of <j. For Sco X-1, GX 17+2, Cyg X— 2 and GX 340+0, the model consists of two COMPTB 
components: one represents pure thermal Comptonization {dotted- dashed line), the other one represents thermal plus bulk Comptonization 
{long-dashed line). For GX 3+1 and GX 354-0 the model is a simple COMPTB with pure thermal Comptonization (solid line). Gaussian 
emission lines, when required by the data, are also shown (dotted line). 



responding to no direct seed photons emission, see Eq. 
[3]) and S=Q (no bulk contribution); hereafter we caU 
this COMPTB as thermal COMPTB. A second COMPTB is 
used instead of the bb and PL components, letting, this 
time, \og{A) and S free to vary during the fit (Table 2). 
The main idea is that the first COMPTB should just re- 
produce the TC spectrum from the outer TL and pre- 
sumably most of a Comptonized disk component, while 
the second one takes into account the TC and BC of the 
NS surface seed photons, which produce the observed 
PL-like hard X-ray tail. For the low/hard state spec- 
tra, GX 354-0, and the high/soft state spectra, GX 3-1-1, 
that are characterized by pure thermal Comptonization, 
the COMPTT model was replaced by a single thermal 
COMPTB (Table 2) with no need for the second, bulk- 
related, COMPTB component. For all of the sources, the 
optical depth r of the TC region is later estimated using 
the COMPTB best-fit values a and kT^ and equations (17), 



(24) in TL95 and then compared to the value reported 
from COMPTT (see Table [1]) . 
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TABLE 1 

Best-fit parameters to the analyzed sources using comptt as thermal Comptonization model. For GX 17+2, Cyg X-2 and 

GX 340+0 additional bb and PL components are required, while for SCO X-1 only a PL IS REQUIRED TO DESCRIBE THE 
high-energy part of the SPECTRUM. FOR GX 3 + 1 AND GX 354-0 THE CONTINUUM MODEL IS A SIMPLE COMPTT. EMISSION LINES, 
WHERE OBSERVED, ARE MODELED WITH A SIMPLE GAUSSIAN MODEL. PARAMETERS WITHIN SQUARE BRACKETS ARE KEPT FROZEN IN THE FIT. 
PHOTOELECTRIC ABSORPTION IS COMPUTED USING THE WABS MODEL IN XSPEC. ERRORS ARE COMPUTED AT 90% CONFIDENCE LEVEL FOR A 

SINGLE PARAMETER. 



Parameter 


Sco X-1 


GX 17+2 


Cyg X-2 


GX 340+0 


GX 3+1 


GX 354-0 




[0.18] 


^■^'-'-0.26 




7 92+012 

' -^^-0.20 


[1.7] 


[3.0] 


BB 


kT^i^ (keV) 
r[2 (km) 




-1 =7+0.09 
^■^'-0.08 


1 10+0.03 
-'-lS-0.04 


1 40+°°* 
^^■^"-0.05 

"•0-0.5 






COMPTT 


fcTi,'*' (keV) 
fcTc (keV) 

T 


J^-U»_0.07 
r, 7=+0.03 
^•'•^-0.03 

4 9+0.2 


n 07+0.07 
■J-'^' -0.06 
= n + 0.2 

a.y_(j 2 


Q 17+0.03 

'-'•^'-0.03 
2 70+0.04 

^- '''-o.oe 
D.y_Q 2 


n 4Q+0.09 
"■^■^-0.33 
2 on+0.04 
^■»0_o 04 

9 2+''-* 
^•^-1.0 


0-62l°;°? 

2 54+0-10 

= q+0.27 
"■•-'-0.30 


1 10+0. 04 
1-18-0.04 

31.lli3^« 

1 =+0.6 
^■"-0.6 


POWERLAW 


r 

pi 


4 1K+0.08 

= 24+152 


9 7 = + 0.15 
^- '''-0.23 

1 7= + 1.19 
^■'°-1.04 


2 =2+0.21 
^•"^-0.21 

f. 00+0. 30 
"■•^'^-0.17 


4-1+0.05 

•^■^J^-0.05 
"•^-1.8 






GAUSSIAN 


El (keV) 
ai (keV) 

r(/) 

EWi (eV) 




"■'"-0.05 
24+° " 

5 1+1-° 

4S+1 
+o_j4 


f. 70+0.10 
f-'O-0.09 

21+°-- 
'-'■^-'--0.15 


6 7Q+0'07 
"■ '^-0.07 

29+°'l'' 
"■^^-0.14 

0+0.8 
'-'■•^-0.8 


[6.4] 

= 7+104 
-126 




xVdof 


212/167 


155/144 


196/152 


174/137 


50/42 


151/131 



^ Equivalent hydrogen column in units of 1022 cm~2. 
^ Temperature of the blackbody spectrum. 

Computed assuming isotropic emission and given distance for each source (see text). 

Temperature of the Wien seed photon spectrum of COMPTT. 
" Powerlaw normalization in units of photons cm~2 §— 1 keV~i at 1 keV. 
' Total photons in the line in units of 10~^ cm~2 s~l. 



3.1. Sco X-1 

The spectrum of Sco X-1 was obtained from observa- 
tions of the INTEGRAL satellite (Winkler et al. 2003) 
using data from the Joint European X-ray Monitor 
(JEM-X, Lund et al. 2003) and the high-energy imager 
IBIS/ISGRI (Ubertini et al. 2003, Lebrun et al. 2003). 
The ISGRI spectrum was obtained by summation of all 
the INTEGRAL observations of the source in the period 
2003-2005 (see P06), while for JEM-X, because of the 
very high brightness of the source (which may give rise 
to problems in the standard observing mode of the in- 
strument) we considered a single pointing (Tcxp ''^ 1800 
sec). We thus caution the reader that the high-energy 
source spectrum (> 20keV) shows, in fact, the average 
behaviour of the source. 

The inclusion of the JEM-X data at lower energies is 
very important in order to get a wide energy coverage, 
which is of key importance to constrain the model param- 
eters. Using these two instruments we can indeed ana- 
lyze the source spectrum in the energy range 4-150 keV. 
The lack of data below 4keV however, does not allow 



XSPEC to determine the photoelectric interstellar absorp- 
tion A^H along the source direction, so we fixed its value 
to 0.19x10^^ cm~^, as found by the Sohd-State Spec- 
trometer on-board the Einstein Observatory (Christian 
& Swank 1997). A simple COMPTT plus PL model for- 
mally would not be completely acceptable from a sta- 
tistical point of view (x^/dof=212/167), even though no 
systematic deviation are observed in the residuals be- 
tween the data and the model. Moreover, this result is 
obtained adding a systematic error of 1% to both ISGRI 
and JEM-X spectra, which is likely below the current 
uncertainty on the instrument calibrations. If system- 
atics are increased to 2%, we obtain x^/dof=189/167. 
The best-fit parameters of the COMPTT plus PL model 
are reported in Table [T] 

We subsequently use a two-COMPTB model: in this 
case, we deal with two seed photon temperatures, which 
can not be determined simultaneously, because the JEM- 
X low-energy threshold 3 keV) is well above the peak 
of the softer seed photon temperature usually observed 
in LMXBs (-0.5 keV, see Barret 2001). 

We thus fix kTs related to the thermal COMPTB to 
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TABLE 2 

Best-fit parameters of the comptb model to the analyzed sources. When fitting the source spectra with two comptb 
components, their electron temperatures kt^^ (thermal plus bulk component) and fct* (thermal component) are kept 
equal each other during the fit. emission lines, where observed, are fitted with a simple gaussian model. parameters 
within square brackets are kept frozen in the fit. parameters within square brackets are kept frozen in the fit. 
Photoelectric absorption is computed using the wabs model in xspec. Errors are computed at 90% confidence level for a 

single parameter. 



Parameter 


Sco X-1 


GX 17+2 


Cyg X-2 




GX 340+0 


GX 3+1 


GX 354-0 




[0.15] 










[1.7] 


[3.0] 








COMPTB (thermal: log(A): 


=8, <5= 


=0) 






fcTi*"^ (keV) 
fcT* (keV) 
a 

caf('=) 


[0.4] 

2 66+° °'^ 

"■^°-0.04, 
2.1 


q qO+0.07 

0.00_g Qg 

n Q'5+0-04 

U.aJ_g Q4 

2.1 


16+"-°'' 

r, 7^+0.05 
^■'4-0.04 

3.4 




Z.aO_g Qg 

n sfi+0-04 

U-OD-0.03 

2.5 


61+° i° 

1 26+°ll 
1.5 


261* 

1 fiq+O-08 
^^■"^-0.08 

1.7 








COMPTB (thermal plus 


bulk) 








fcTi''' (keV) 

log(A) 

a 

6 


1-^1-0.10 
-n 47+0-52 
o (;o + 0.07 

>63 

i=kn] 


1 45+0.30 
l-43_0.26 

0.28 (>-1.00) 
2.37 (< 2.73) 
89 (>15) 

i=kn] 


0.26 (>-0.14) 
[2.5] 

60 (>27) 

i=kn] 




-I qO+0.18 

_0 44+0-54 

"^-0.55 
52 (>23) 

i=kn] 






GAUSSIAN 


El (keV) 
o-i (keV) 

Ad) 

^1 

EWi (eV) 




6 69+° °5 
D-Da_0.05 

24+°-l° 
"-'•''4_o.09 

5 5+1-3 

47+3 
4'-14 


2 5+1° 
39118 




6 73+0-07 
°- ''^-0.10 

n S4+0-26 

4 3+2.19 
^■•^-0.90 


[6.4] 

0.95(> 0.51) 
286l?^g 






3.65 


1.53 


0.96 




1.42 


0.19 


0.06 


T f"^' /T 
-^comptb, th/^tot 


0.79 


0.86 


0.74 




0.90 








214/165 


158/143 


189/152 




182/136 


48/42 


155/131 



" Equivalent hydrogen column in units of 10-^^ cm~2, 

^ Temperature of the COMPTB blackbody seed photon spectrum (same as fcT^b in Tab. [TJ. 

Compton amplification factor, see text for definition. 

Total photons in the line in units of 10~3 cm~2 s~^. 
" Estimated in the energy range 0.1-200 keV. 



0.4 keV, allowing instead to vary the seed photon tem- 
perature kTs suggested to be related to the NS surface. 
The results of the fit with a two-COMPTB model are re- 
ported in Table [2l from the best-fit values a and kT^ of 
the thermal comptb component and equations (17) and 
(24) in TL95, we estimate r=5.6, which is not far away 
from that directly find by COMPTT ('^ 5, see Table [T]). 
As it concerns the second COMPTB component, aimed 
to describe the hard X-ray tail, we find that the bulk 
parameter d would be pushed up to unreasonably high 
values by xSPEC, with no apparent minimum in the 
space. It is, in fact, just possible to put a lower limit on 
S, as reported in Table O whereas good constrains are 
instead obtained for the seed photon temperature kTs 
and the energy index a. In Figure [3] we report the de- 
convolved unabsorbed source spectrum and best-fit two- 
COMPTB model. For estimating the source 0. 1-200 keV 
luminosity we assumed a distance of 2.8 kpc (Bradshaw 
et al. 1999). 



3.2. GX 17+2 

For GX 17+2 we use BeppoSAX (Boella et al. 1997a) 
observations performed on April 3, 1997, when the source 
was in the left part of the horizontal branch (HB, see Fig. 
2 in F07) and a hard X-ray tail was present in the source 
spectrum. The data set is the same previously used in 
the analysis by F05 and F07. 

Following Di Salvo et al. (2000b) and F05, we fit the 
0.4-120 keV source spectrum with a photo-electrically 
absorbed bb plus TC (comptt) component, plus PL to 
fit the hard X-ray tail, plus a GAUSSIAN emission line 
around 6.7 keV. The best-fit parameters of the model are 
reported in Table[TJ Assuming a distance from the source 
of 7.5 kpc (Penninx et al. 1988) and using the best-fit bb 
temperature and fiux, we estimate a bb radius, ^ 5 km. 
While keeping in mind all the limitations in making such 
an estimate (i.e. color to effective temperature correc- 
tion, possible anisotropy in the emission), the order of 
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magnitude of the bb radius supports the idea that this 
component very Ukely originates close to the NS surface. 
The PL photon index T is the same aheady reported in 
Di Salvo et al. (2000b) and F05. 

In the next step, we have good constrains on the best- 
fit parameters of the pure thermal COMPTB component 
using the two-COMPTB model: they are reported in Ta- 
ble [2] Again, from the thermal a and /cTo-values, we 
estimate r=5.2, close to r ~ 6 found by COMPTT. The 
situation is a bit more critical in the thermal plus bulk 
related COMPTB component. In particular, for the illu- 
minating factor log(j4) and the bulk parameter 6, the 
lower limits (at 90% confidence level) can be only put by 
XSPEC, while, on the contrary, we obtaine only the up- 
per limit on the a-value. These statistical uncertainties 
are mainly related to the small amount of high-energy 
(> 30keV) points, as it can be seen in Figure O where 
we report the deconvolved source spectrum with the two- 
COMPTB model. 

3.3. Cyg X-2 

For Cyg X~2 we use a BeppoSAX data set which has 
been already analyzed by DS02. These data are related 
to a source observation performed on 1996 July 23. Dur- 
ing this observation, the source traced the full HB across 
its colour-intensity diagram. The HB was divided into 
two parts (upper and lower HB, respectively, see Fig. 2 
of DS02) for which separated spectral analysis was car- 
ried out; in both cases a hard X-ray tail was detected. 

In this Paper we report results from the upper HB. 
DS02 fitted the 0.1-lOOkeV continuum with a niuhi- 
colour disk blackbody (diskbb in XSPEC, Mitsuda et al. 
1984) plus COMPTT plus PL plus two GAUSSIAN emis- 
sion lines at ~ IkeV and ~ 6.7 keV, respectively. We 
note however that DS02 used for the Low-Energy Con- 
centrator Spectrometer (LECS, Parmar et al. 1997) the 
standard on-line available response matrix. The latter 
may give rise to spurious features when the instrumental 
count rate exceeds ~ 50 counts s^^. For this Cyg X-2 ob- 
servation, the LECS count rate is ~ 60 counts s"'^, thus 
requiring the production of an observation-related re- 
sponse function, which was generated using the LEMAT 
package (vS.O.l). 

With the newly produced matrix we establish, in fact, 
that there is no more strong evidence of the 1 keV emis- 
sion line. On the other hand, an excess around 2.6 keV 
is still found in the residuals of both the LECS and 
the Medium-Energy Concentrator Spectrometer (MFCS, 
Boella et al. 1997b). We fit this excess with a Gaussian 
emission line, even though it is not clear whether its ori- 
gin may be attributed to the source or arise from some 
instrumental effect. The results of the fit with a bb plus 
COMPTT plus PL model are reported in Table [T] The PL 
index F results steeper than that reported in DSG2. This 
is mainly due to the different response matrix used at 
low energies (< 4keV), namely in the region where the 
interstellar absorption A^h and PL normalization must be 
simultaneously determined by xspec. We also estimate 
the BB radius, finding ~ 8 km, for a source distance of 
8kpc. 

Replacing the above model with the two-COMPTB one, 
similarly to Sco X-1 and GX 17-f2, we find very good 
constraints on the best-fit parameters related to the ther- 
mal COMPTB (see Table [2|). As in the previous sources, 



the optical depths inferred from best-fit a and fcTc- 
values, and estimated by COMPTT are very close each 
other, (r ~ 6.3 and r ~ 6.9, respectively). 

As it concerns the thermal plus bulk parameters, unlike 
Sco X-1 and GX 17+2, where, even though poorly con- 
strained, they can be left free in the fit, there is no way to 
perform serious analysis leaving all the parameters free. 
For Cyg X-2 we thus performed spectral fitting in sev- 
eral steps with a fixed at each step in the range from 0.5 
to 3.5, finding that the minimum in the parameter 
space is obtained for a=2.5. Fixing a, just allows us to 
put lower limit (90% confidence level) for 5 and \og{A). 
In Figure [3] we report the deconvolved source spectrum 
with the two-COMPTB models. 

3.4. GX 340+0 

The spectrum of GX 340-1-0 refers to a pointed 
BeppoSAX observation of the source performed between 
2001 August 9 and 10. A more detailed description of 
the data reduction can be found in L04. The results 
with the BB plus COMPTT plus PL model, with an addi- 
tional GAUSSIAN emission line, are reported in Table [TJ 
It is worth pointing out that the COMPTT model with 
slab geometry gives a slightly worse result with respect 
to that with spherical geometry (for which we obtain 
X^/dof=:160/137). We also want to emphasize the dif- 
ferent values of the direct BB and COMPTT seed pho- 
tons temperature (fcTbb and kT^, respectively) reported 
in Table [T] in comparison with the results of L04. As 
discussed in F05, the usually adopted model bb plus 
COMPTT for description of the persistent continuum spec- 
tra of LMXBs, may give rise to a dichotomy in the so- 
lutions, given that two solutions (statistically indistin- 
guishable) are possible, one with fcTbb > kT^^ and the 
other with fcTbb < A:T„. In the former case, the BB-like 
emission is suggested to mainly come from the NS sur- 
face, while in the latter one it is presumably related to 
the accretion disk. These different interpretations are 
also suggested by the inferred BB-radii which are of or- 
der of the NS radius in the former case or significantly 
greater 78 km in L04) in the latter one. On the basis 
of the accretion scenario proposed by F07, we prefer the 
first solution (see Table [T]). 

The results with the two-COMPTB model are instead re- 
ported in Tabled! In this case, the general considerations 
on the other three sources can be extended to GX 340-t-O: 
the constraints on the best-fit parameters of the thermal 
COMPTB component are very good, while they are worse 
for the thermal plus bulk related one, except for the seed 
photon temperature. Noticeably, the inferred r from the 
best-fit (thermal) a and kT^ is lower (r ~ 6) than the 
one directly estimated with COMPTT (r ~ 9, see Discus- 
sion). For the estimation of the bolometric (0. 1-200 keV) 
source luminosity we assumed a distance of lOkpc (L04). 

3.5. GX 3+1 

We use the data from the RXTE (Bradt et al. 1993) 
public archive related to a source observation performed 
on August 21 2004. A joint spectrum from the second 
units of the Proportional Counter Array (PGA, Jahoda 
et al. 2006) and High-Energy X-ray Timing Experi- 
ment (HEXTE, Rothschild et al. 1998) on-board the 
spacecraft was obtained, following the standard proce- 
dure for data reduction. At high energies the source 
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was detected by HEXTE only up to 30 keV. Similarly to 
See X-1 INTEGRAL spectra, the low-energy threshold 
of 3keV does not allow us to constrain the photoelec- 
tric interstellar absorption A'h in the source direction, so 
we fixed its value to 1.7x10^^ cm^^ (Christian & Swank 
1997). 

We find that a simple photoelectrically-absorbed 
COMPTT model is good enough to describe the source 
continuum spectrum. However, a clear excess in the re- 
gion 6-7 keV reveals the presence of a iron emission line 
which is fitted with a simple GAUSSIAN. The poor instru- 
mental resolution in this energy region, and the uncer- 
tainty calibration of the instrument because of the diffuse 
Galactic ridge emission, do not allow us to simultane- 
ously constrain all the line parameters. We thus fix the 
line energy centroid at 6.4 keV during the fit. Because 
of these uncertainties, both the line broadening ai and 
equivalent width values, must be considered with some 
caution. The best-fit parameters of the COMPTT plus 
GAUSSIAN model are reported in Table [TJ 

As in the case of Sco X-1 fsee 13. ip . the absence of a 
direct BB-like component in the spectrum is very likely 
due to a combination of the low-energy PCA thresh- 
old (3keV) and data quality. Subsequently, replacing 
COMPTT with thermal COMPTB, we find the same statis- 
tical result, with best-fit values reported in Table [51 The 
inferred optical depth of the TC region, obtained from 
best-fit kTc and a values, is r ~ 4.9, very close to that 
found by COMPTT (t ~ 5.3, see Table [T]). The bolomet- 
ric 0. 1-200 keV source luminosity is estimated assuming 
a distance of 5 kpc. We note that this value is certainly 
underestimated, as it includes the region below 3 keV, 
not covered by PCA, where a BB-like component has 
been already observed from the source (Oosterbroek et 
al. 2001). 

3.6. GX 354-0 

The spectral analysis on GX 354-0 is carried-out using 
data obtained by F06 with JEM-X and ISGRI onboard 
the INTEGRAL satellite. Among the nine spectra ex- 
tracted by FOG as a function of the source position in the 
HID, we consider the hardest one (see Fig. 3 and Table 
2 in F06). Unlike FOB however, we assume a spherical 
geometry for the plasma in COMPTT, given that it seems 
to be more representative for the hard state accretion 
scenario. A simple photoelectrically-absorbed COMPTT 
model is good enough to fit the source spectrum and no 
soft X-ray emission is seen in the spectrum. 

Replacing COMPTT with a thermal COMPTB, we ob- 
tain the same results of COMPTT, both in terms of statis- 
tics and as it concerns the agreement between the seed 
photon and electron temperatures (see Tabs. [T]and[2]). 
Using equations (17) and (24) of TL95 for the spheri- 
cal case, from the best-fit kT^ and a-values we derive 
T ~ 1.5, which perfectly matches the best- fit value di- 
rectly obtained from COMPTT. The source bolometric 
0. 1-200 keV luminosity was computed assuming a source 
distance of 5 kpc (Di Salvo et al. 2000a; Galloway et al. 
2003). Figure [3] shows the deconvolved source spectrum 
and superposed best-fit model. 

4. DISCUSSION 

We present a new model for the X-ray spectral fitting 
package XSPEC, named COMPTB. The principal aim of 



this model is to provide a more physical description of 
the transient PL-like hard X-ray (> 30 keV) tails observed 
in LMXBs, in the framework of the bulk motion Comp- 
tonization theory. We recall however that COMPTB is 
actually not just specific to bulk motion. Indeed, setting 
the bulk parameter 6 to zero, it reduces to a generalized 
thermal Comptonization model, in the same fashion as 
the widely used models COMPST or COMPTT. 

4.1. Thermal Comptonization component 

In the six sources that we analyzed, the spectra are 
dominated by a TC component which carries-out most 
of the total energetic budget (see Tableland Fig. [J). 
In the case of GX 34-1 and GX 354-0 this is in fact the 
only component observed in the X-ray spectrum. 

In the GX 3+1 case, this is a result of a low-energy 
(~ 3 keV) threshold effect of PCA, given that a BB-like 
emission has been previously observed with BeppoSAX 
(Oosterbroek et al. 2001). On the other hand, in the case 
of GX 354-0, it is possible that this soft emission was in- 
trinsically absent or very weak in the analyzed spectrum 
of the source given that it was detected when the source 
was in the intermediate state (see F06 and also Di Salvo 
et al. 2000a). F06 attributed the appearance of the soft 
BB-like emission, in the soft state, to the accretion disk. 
Instead, in our accretion scenario (as already reported 
in F05), the BB-like emission is related to the emission 
from the NS star surface. Independently of the emission 
region from which it comes, the fact that this compo- 
nent is not observed in the low/hard state spectrum of 
GX 354-0 is consistent with an accretion geometry where 
the TC corona is significantly extended and being almost 
spherical. This spherical configuration should completely 
cover the central object, intercepting all the seed photons 
coming from the NS surface. 

For all of the studied sources, we compute the Comp- 
ton amplification factor (CAE) defined as the ratio of the 
Comptonized to seed photon energy fiuxes. Noticeably, 
the Compton amplification factor (CAE) gets the highest 
value in Cyg X-2 (3.5), despite the fact that its TC com- 
ponent provides the lowest fractional contribution to the 
total spectrum among the six studied sources (see Table 
[2]). This is ultimately related to the fact that the seed 
photon energy is softer (fcTg ^ 0.2) in Cyg X-2 than in 
the remaining five sources (see Chakrabarti & Titarchuk 
1995, for the CAE dependence on the seed photon en- 
ergy). 

It is important to emphasize that the energy index 
a directly provides the efficiency of the Comptoniza- 
tion process (see e.g. Bradshaw et al. 2007, here- 
after BTK07), no matter which is the plasma geometry 
(i.e. slab or sphere). For a given geometry (taking into 
account hydrodynamical considerations), one can infer 
Thomson optical depth r using the best-fit a and kT^ 
values provided by COMPTB and using equations (17) and 
(24) in TL95. 

We find that the derived values are in very good 
agreement with those obtained from COMPTT, except for 
GX 340-1-0, where r provided from COMPTT (~ 9.2) is 
higher than that estimated using COMPTB 5.8). This 
difference likely comes from the slightly different descrip- 
tion of the Comptonization by these two models. 

4.2. Combined thermal-bulk Comptonization component 
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The bulk motion Comptonization process as a possi- 
ble origin of hard X-ray tails in LMXB systems hosting 
either a NS or a BH has been suggested from both ob- 
servational (Shrader & Titarchuk 1998; Borozdin et al. 
1999; P06, F07) and theoretical (TMK96, TMK97, Lau- 
rent & Titarchuk 1999, 2001; TF08) works. A model 
which can be used to describe these emerging spectra 
is already present in the xspec package (bmc, TMK97), 
but as previously discussed in Section [U it presents some 
limitations, which are the lack of the electron recoil term 
(cutoff) in the Green's function and the parameteriza- 
tion of the high-energy Comptonized component in terms 
of just the spectral index a (no S bulk parameter). In 
fact, index a has only information about the efficiency of 
Comptonization, no matter whether it is purely thermal 
or due to combined effects of thermal and bulk Comp- 
tonization (see TMK97 and BTK07). 

We try to overcome these limitations, including the 
full Green's function expression (Eq. [3]) in our model 
and using numerical convolution of the Green's function 
with the seed photon spectrum (Eq. [5]). 

One of the advantages of this numerical approach is 
that the condition kTs <C Epi (seed photons energy much 
less than plasma energy and Green's function approx- 
imated by a broken PL, see ST80, TMK97) is not re- 
quired. This model modification is, in fact, important 
in the spectra of high- luminosity NS LMXBs, where the 
seed photon and plasma energy (electron temperature) 
differ just for a factor of a few (see Tables [T] and [21) . 

We apply our combined thermal plus bulk Comp- 
tonization model {S > 0) to the four LMXBs of our 
sample which shows a PL-like hard X-ray tail, namely in 
Sco X-1, GX 17+2, Cyg X-2 and GX 340-1-0 spectra. In 
particular, we use a two-COMPTB model for these sources 
in which the first COMPTB takes into account the pure 
TC part of the spectrum (see previous subsection) while 
the second COMPTB describes both the BB-like emission 
observed at low energies and the PL-like emission above 
SOkeV. We assumed the same temperature kT^ for the 
thermal and thermal plus bulk Comptonization region. 
The correctness of this assumption depends on how the 
two regions are physically and geometrically coupled and 
whether an actual temperature gradiend does exist on 
the region where Comptonization occurs. 

In the accretion geometry proposed by F07, bulk 
Comptonization occurs in the innermost part of the tran- 
sition layer region, while TC is dominant in the outer TL 
and presumably in some extended region located above 
the accretion disk. The latter statement is suggested by 
the lack of a direct soft (~ 0.5 keV) BB-like X-ray emis- 
sion in the spectra, usually attributed to the accretion 
disk. It is likely that most of the disk emission is up- 
scattered in the TL or somewhat embedded in the total 
spectrum, without possibility to split it as a single spec- 
tral component. We find some indication of contribution 
of disk soft seed photons to the TC component using 
the best-fit values of the BB temperature of these pho- 
tons. Note this temperature is significantly lower than 
that of the BB-like seed photon component presumably 
coming from NS surface (see Tabs. [T] and [2] for COMPTT 
and COMPTB, respectively). Particularly for Cyg X-2 we 
obtain a value < 0.2 keV, in contrast to higher values 
of that found by DS02 or by Done & Gierhnski (2003). 
Higher disk temperatures are also found by Gilfanov et 



al. (2003) and Revnivtsev & Gilfanov (2006). 

However we should point out that all these authors 
(but DS02) used PCA/RXTE data for which the low- 
energy band is higher than 3 keV, whereas our best-fit 
seed photon temperatures are based on the BeppoSAX 
spectra in the broad energy band from 0.4 to 120 keV. 
The PCA/RXTE with its low-energy threshold above 3 
keV cannot in principle determine any trace of the soft 
component whose temperature is lower than 1 keV given 
that its maximum emission (< 3 keV) would still fall 
outside the PCA energy band. On the other hand the 
disk photon temperatures of 0.8-1.4 keV found by DS02 
is inferred using a response matrix for the low-energy 
instrument (LEGS) which is not correct for such bright 
sources and thus (as also testified by strong decrease of 
the 1 keV emission line in our data) must be considered 
very cautiously. 

The uncertainties of the bulk-related parameters of 
COMPTB (see Table [2]) do not allow us to check the ac- 
tual differences among them as it concerns their inner- 
most bulk-dominated region given the presence of a few 
points only at high-energies (> 30 keV) in the spectra 
of the four intermediate state studied sources (see four 
upper panels in Fig. [3]). 

However we should note that in all cases 5, through 
which the relative importance of BC to TC is 
parametrized, is of order of tens, which values are ex- 
pected in the case of efficient BC, when the hard PL-like 
emission should be seen. The fact that we can not con- 
strain the upper value (at 90% confidence level) of S is 
not surprising because of the measurement of the high- 
energy cutoff (not observed in the analyzed spectra) is 
required. 

Indeed, the most important effect of the bulk term in 
the Green's function is to push forward the rollover en- 
ergy, with respect to the pure TC case: higher (5-values 
lead to higher energies at which the spectrum exponen- 
tially falls down. This can be clearly seen in Figure [TJ 
On the other hand, it is possible to provide only a lower 
limit of S if the high-energy behaviour of the available 
data goes as an unbroken PL. This is exactly what hap- 
pens for the INTEGRAL/miS spectrum of Sco X-1 (see 
Fig. ED. 

In the case of GX 17+2, Cyg X-2 and GX 340-hO, the 
thermal plus bulk COMPTB component predicts a rollover 
of the spectrum in the energy range 100-200 keV (see Fig. 
[3]). However, it is not possible with the present data to 
really confirm this prediction, even the fact that xspec 
finds a minimum of S in the parameter space could 
not be indicative of that. 

In the thermal plus bulk COMPTB component, the 
Green's function energy index a is, for all sources, higher 
than that obtained for pure TC, with the highest dif- 
ference observed in Sco X-1. In Figure [U we plot the 
Green's function for monochromatic input line at energy 
a;o=0.3 (where xq = Eo/kTf.), for two different values 
of a (1 and 2.5, respectively) considering pure thermal 
((5=0) and bulk-dominated (J=30) cases. The chosen 
value of Xq is, in fact, representative of the NS local en- 
vironment condition under this study, as kTg ~ 1 keV and 
fcTo ~ 3keV (see Tables [U and d]). Looking at Figure [J 
one can note that high (5- values push forwards the high- 
energy cutoff, as already mentioned above. It is evident 
that the higher a-values related to the BC region with 
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respect to those related to TC region (see Table [T]), are 
indicative of the fact that in the innermost part, for all 
four sources, the TC efficiency is reduced but the BC ef- 
fect dominates over the TC one given that the power-law 
hard tail is detected. 

4.3. Limitations and applicability of the model 

When one uses COMPTB to fit the X-ray spectra of NS 
LMXBs, there are some issues which must be kept in 
mind. 

Let us first consider the case of pure TC spectra {6—0 
in the Green's function, see TMK97 and Eq. [3] here). 
The Green's function of the model is derived as a solu- 
tion of the diffusion Comptonization equation, which is 
in turn obtained as an approximation of the Boltzmann 
kinetic equation in the diffusion regime. This means that 
the average photon energy exchange per scattering is 
small {AE/E ^ 1) and photons suffer many scatterings 
when traveling across the plasma before escaping. 

This condition holds when the plasma temperature is 
about 20-30 keV and less but its optical depth is high, 
such as observed in the persistent spectra of Z sources 
and high-luminosities atoll sources (see Barret 2001 for a 
review) . In general most of the NS LMXBs obey this con- 
dition, as the strong radiation ffux coming from the cen- 
tral object acts as a thermostat in controlling, through 
Compton cooling, the plasma temperature (Barret et al. 
2000). 

In the case of BHCs, the situation is quite different and 
typical hard states for this class of sources are character- 
ized by plasma temperatures of 60 keV (see e.g. Remil- 
lard & McClintock 2006) and optical depth a few. In this 
case, the diffusion approximation should be used with 
some care (the photon energy exchange for scattering is 
not so small but the number of scatterings before pho- 
tons escape is relatively small), and a different treatment 
of the Compton scattering should be addressed: models 
such as COMPPS (Poutanen & Svensson 1996), which use 
iterative scattering numerical procedures are thus most 
suitable to be used for fitting BHCs hard state spectra. 
This is one of the reasons why COMPTB (but also COMPST 
or COMPTt) should be apphed preferably to NS LMXBs 
spectra with plasma temperatures below 20-30 keV. In 
other cases, BHs or NS in their hard/PL state (see Sec- 
tion 3), the model can just fit the shape of the spectrum: 
while this may result satisfactory in terms of value, 
the output parameters should be taken cautiously. 

Let us now concentrate on the case where BC domi- 
nates over TC {6 ^1). The analytical solution of the 
Comptonization equation ([T]), for the Green's function, 
does exist only neglecting the term, (wb/c)^/(30) in 
equation ([T]), where 8 = kTc/rrioC^. Whether this as- 
sumption is correct enough, it depends on the local en- 
vironment conditions close to the central object. In the 
case of systems hosting a NS, the presence of a firm sur- 
face plays a key role in determining the global inward 
motion of the accreting material. At high accretion rates, 
the NS surface on one side is the source of a strong radia- 
tion pressure due to kinetic energy release of the matter, 
but on the other hand its mirror-like behaviour (reflect- 
ing inner boundary condition) is crucial in determining, 
at any radius, the pressure gradient due to local gravita- 
tional energy release of the in-falling material. 

In many cases, these concurring effects may efficiently 



stop (or strongly decrease) the in-falling matter velocity 
fb, thus suppressing the bulk Comptonization effect (and 
associated hard X-ray tail). Numerical simulations by 
TF08 show that under this conditions Ub is ^ 0.2 c, and 
typical in-falling velocities close to the NS surface are ~ 
0.1 c. In this case, for an electron temperature kTc=3 
keV, it results (ub/c) 7(39)^0.6. 

In the case of BHCs, the different inner boundary con- 
dition at the event horizon (e.g, Titarchuk & Fiorito 
2004) strongly modifies the radiation pressure behaviour 
with respect to the NS case, and may lead the inward 
bulk velocity Vh to get values of 0.3-0.4 c in proximity 
of the last stable marginal orbit, even at accretion rates 
higher than Eddington limit. Given that the typical tem- 
peratures for these states are similar to those of the NS 
intermediate and very soft states (e.g., Titarchuk & Fior- 
ito 2004), namely ~ 3-5 keV, we derive that ivb/cf/{3Q) 
may reach a few. The full kinetic treatment should be 
used to determine the Green's function for this case (see 
Titarchuk & Zannias 1998, Laurent & Titarchuk 1999) 
because the value of Ubii0.3c and the photon energy 
change is not small in the bulk flow into BH. 

The numerical solution of the Comptonization equa- 
tion including the second order bulk term shows that it 
produces two net effects on the the emerging spectrum: it 
pushes forward the high-energy cutoff and makes higher 
the Comptonized spectrum normalization (see Fig. 1 in 
TMK97) . In terms of our model parameterization, given 
that the high-energy cutoff is not observed in the data, 
we would expect that this higher Compton normaliza- 
tion, actually due to (wb/c)^, would be somewhat com- 
pensated by a higher illuminating factor \og{A) value. 
The model can thus always fit the shape of the spec- 
trum, whereas its best-fit parameters should be taken 
with some precaution as in the case of pure thermal, 
high-fcTc, TC spectra (see above). 

A second issue to be pointed out is how one should de- 
fine 6. The analytic form reported in equation ^ is de- 
rived for a plasma subjected to a pure free-fall radiation- 
corrected velocity profile, where only gravity and radia- 
tion pressure from the central object are the concurring 
processes. In particular, if the compact object is a BH, 
only gravity is present. 

The condition where only gravity and radiation pres- 
sure are at work holds, in fact, in the case of a optically 
thin medium, where the radiation fiux (from the NS) 
follows the law. However, it is known that high- 
luminosity LMXBs systems which show hard X-ray tails 
(the six known Z sources and, recently, also GX 13-f 1, 
see P06) are characterized by high M values (mainly in- 
ferred from their luminosity), and consequently optically 
thick environments. In this case, additional (and even 
dominant) effects due to radiation pressure from local 
gravitational and kinetic energy releases, viscous energy 
transportation, gas and magnetic pressure, cannot be ne- 
glected. The numerical simulations performed by TF08 
show that under these conditions the radial velocity be- 
haviour of the in-falling material deviates from a pure 
free-fall law across the whole TL. Nevertheless, in the 
innermost part of the TL itself, gravity and radiation 
pressure actually dominate and the velocity profile is al- 
most free-fall radiation-corrected like and d definition as 
defined in equation ([2]) can be almost safely used. 
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To strengthen this statement, we tried to make a in- 
dependent (5-estimation as it would be derived in the 
opticaUy thin case (see Appendix). It is worth not- 
ing that in equation (|A3|) . the only parameter which 
cannot be directly obtained from the observations is 
the dimensionless accretion rate m. We find the S- 
values 155/m, 145/rri, 178/m and 169/rh for Sco X-1, 
GX 17-^2, Cyg X-2 and GX 340-^0, respectively. While 
keeping in mind the issues discussed above, the similarity 
of the (5-values, as order of magnitude, between these val- 
ues and those reported in Table [2] is interesting. Higher 
values of m > 1 would even reduce this difference. 

4.4. Similarities and differences between COMPTB and 
BW07 model 

As already pointed out in the Introduction, BW07 de- 
veloped a TC and BC model for accreting X-ray pulsars. 
The radiative transfer form reported in equation (34) of 
BW07 is sligthly different from that reported in equation 
(14) of TMK97, which we adopted in our Paper (see Eq. 
[T]). The difference mainly arises because of the differ- 
ent assumptions on the velocity profiles in TMK97 and 
BW07. 

A free-fall velocity profile v{r) oc r^^/^ is used by 
TMK97 which leads to r(r) cx r^^/^ and v{t) cx T using 
the continuity equation (see TMK97 and also TMK96). 
This result is obtained keeping in mind that the surface 
through which matter fiows scales as r^. 

Instead, in the model of BW07, plasma motion chan- 
neled by the strong magnetic field is modelized as a flow 
passing through a cylindric column whose area does not 
change with the distance r from the NS polar cap (see 
their Eq. [19]). Using this particular form of the continu- 
ity equation BW07 show that r(r) c>c r^/^. In order to an- 
alytically treat (by mean of variable separation) the prob- 
lem, BW07 assume that v{t) cx r and thus v{t) oc r^^^. 
This direct, instead of inverse, proportionality between 
w(r), r(r) and r would obviously modify the form of the 
spatial term in the Comptonization equation when writ- 
ing it as a function of the optical depth r instead of the 
distance r. 

Moreover we note that in BW07 optical depth r is de- 
fined as such that dT=Ne(J\\dr (where is the Thomson 
cross-section for photons propagating in parallel direc- 
tion to the magnetic field lines and r is a radial length 
coordinate), while TMK97 use an effective optical depth, 
defined such that Toff (rti.ap)=l, where rti-ap is the pho- 
ton trapping radius. Looking at equations (1) and (5) 
in TMK97 is evident that with such a parametrization 

Toff (r) oc v{Tes) oc t^"^^. 

It is worth noting that also 13W07 neglect the second 
order velocity (ub/c)^ bulk term, which is essential in or- 
der to separate Comptonization equation in space and 
energy. Another important difference between the two 
models is that BW07 study the BC effect in the presence 
of a strong magnetic field {B ^ 10^^ G) which modifies 
the electron cross section in the photon direction parallel 
and perpendicular to the magnetic field lines. However 
it is not our case, as the class of sources we consider 
are characterized by B-fields of order B< 10^ G and less 
(see Titarchuk et al. 2001, for details of B-field deter- 
mination in LMXBs). For such low magnetic fields the 
B-field modification of the electron cross-section can be 



neglected. 

However it is remarkable that for both models the 
Green's function is a broken PL in the energy band where 
the recoil effect can be neglected (compare Fig. |3]in our 
Paper and Fig. 5 in BW07). Also it is common for these 
two models that the high-energy PL is followed by ex- 
ponential turnover when the recoil effect is taken into 
account. It implies that the shape of the emergent spec- 
trum is similar and it can be also fitted by the xspec 
models COMPTT and BMC. The main difference is how 
the PL index of the Green's function and high-energy cut- 
off are related to the physical parameters of the particu- 
lar model. We can consider our models to be somewhat 
complementary as, starting from the same physical pro- 
cess but in very different enviroments, they are supposed 
to describe emergent spectra for two classes of NS binary 
systems. 

5. CONCLUSIONS 

We developed a new model for the X-ray spectral fit- 
ting package xspec (comptb) which can be considered a 
general Comptonization model and more specifically the 
extension of the previously developed Bulk Comptoniza- 
tion Model (bmc), already present in xspec. Similarly 
to BMC, our model consists of two components: one rep- 
resents the BB-like seed photons which escape the plasma 
cloud without appreciable energy exchange, the other one 
gives the effectively Comptonized spectrum, obtained as 
the convolution of the system Green's function with the 
seed input photon spectrum. 

Using spectra of the BeppoSAX, INTEGRAL and 
RXTE satellites, we tested the model on a sample of six 
NS LMXBs which experience different spectral states: 
intermediate state (Sco X-1, GX 17+2, Cyg X-2 and 
GX 340-1-0), where a PL-like hard X-ray component ap- 
pears over the stable continuum dominated by thermal 
Comptonization, low-fcTo, high-T spectrum, a very soft 
state (GX 3-1-1), similar to the previous one but without 
the PL-like component, and a low/hard state (GX 354-0) 
where only thermal Comptonization from a hot plasma 
(fcTo ~ 30 keV, T ~ 1.5) is observed in the spectrum. 

We find an excellent agreement, as it concerns fitting 
the TC components, between the best-fit parameters of 
COMPTB and those of COMPTT (the most widely used 
model to this scope), which confirms the goodness of the 
code. 

On the other hand, one of the most important goals 
of COMPTB is to try a more physical approach, in the 
framework of the bulk motion Comptonization process, 
in explaining the transient PL-like hard X-ray (> 30keV) 
emission observed in NS LMXBs. This is what we have 
done in analyzing the spectra of the four intermediate 
state sources. It is worth noting that bulk motion Comp- 
tonization is supposed to be at work, with even higher 
efficiency, also in BH systems in their soft states, where 
a unbroken PL extended up to MeV energy is usually 
observed. 

We show that the values of the (5-parameter, which rep- 
resents the importance of bulk with respect to thermal 
Comptonization, obtained from the fit, can be physically 
meaningful and can quantitatively describe the physi- 
cal conditions of the environment in the innermost part 
of the NS systems. For three of the intermediate state 
sources (GX 17+2, Cyg X-2 and GX 340+0) the best- 
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fit model predicts a high-energy cutoff in the spectrum 
in the energy range 100-200 IceV, which is however not 
possible to confirm on the basis of the present available 
data. 

The next generation of high-energy missions, with 
their improved sensitivity and extended energy cover- 
age, will definitively address this issue. On the other 
hand, the theoretical results coming from a full magneto- 
hydrodynamical treatment of the TL region in NS and 
BH systems are very promising. We are planning to per- 
form Monte Carlo simulations where the derived radial 



velocity profiles for NS systems will be of key importance 
in order to predict the position of the high-energy cut- 
off due to BC effect, to be then compared with future 
observations. 
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APPENDIX 

ESTIMATION OF THE (5-PARAMETER FROM OBSERVABLE QUANTITIES 
The (5-parameter for an optically thin accreting plasma is given by 



(Al) 



where £ = F{R) / FEdd{R) is the ratio of the actual to Eddington fluxes at any distance R, m is the accretion rate in 
Eddington units and O is the dimensionless electron temperature O = kTc/rricC^ . In the case of isothermal plasma, 
given that both F(R) and FEdd(R) go as R^^, the parameter 6 is independent of R. Let us define now i^obs as the flux 
measured by the observer at distance D. This is related to the flux at a given distance R by 



(A2) 



Keeping in mind the definition of the Eddington flux FEdd=GA/c/ic/fci?^ (where fc=0.4 cm^ g"^ and /ic is the average 
electron molecular weight, which we assume equal to 1.14 as for solar abundances) and using equation (|A2[) we may 
rewrite equation (jAip as 



S= 1- 



1/2 



1 

mO 



(A3) 



If we relate Fobs to the BB-like flux coming from the NS, then we can estimate S using the best-fit COMPTB parameters, 
by setting to zero the illuminating factor A, or correspondingly, log(A)^l. This would indeed provide just the flux of 
the direct BB-like component of the spectrum (see Eq. HI) . On the other hand the dimensionless parameter 8 can be 
estimated using the best-fit temperature kTg. Thus, assuming a NS mass M=1.4 Mq and if the distance D from the 
source is known, it is possible to estimate 5 as a function of the mass accretion rate m, which is the only quantity not 
directly measurable from the observations. 
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